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We investigate the stability of theories in which Lorentz invariance is spontaneously broken by 
fixed-norm vector "ffither" fields. Models with generic kinetic terms are plagued either by ghosts 
or by tachyons, and are therefore physically unacceptable. There are precisely three kinetic terms 
that are not manifestly unstable: a sigma model {d^iA^)^ , the Maxwell Lagrangian -Fp^-F*"^, and 
a scalar Lagrangian {dy,A^Y . The timelike sigma-model case is well-defined and stable when the 
vector norm is fixed by a constraint; however, when it is determined by minimizing a potential there 
is necessarily a tachyonic ghost, and therefore an instability. In the Maxwell and scalar cases, the 
Hamiltonian is unbounded below, but at the level of perturbation theory there are fewer degrees of 
freedom and the models are stable. However, in these two theories there are obstacles to smooth 
evolution for certain choices of initial data. 



I. INTRODUCTION 



The idea of spontaneous violation of Lorentz invariance through tensor fields with non- vanishing expecta- 
tion values has garnered substantial attention in recent years jTJ[51|21lll|SJISll71IHll21 EHl EH H] • Hypothetical 
interactions between Standard Model fields and Lorentz-violating (LV) tensor fields are tightly constrained 
by a wide variety of experimental probes, in some cases leading to limits at or above the Planck scale 

[lliaiiallllillsilillllT]. 

If these constraints are to be taken seriously, it is necessary to have a sensible theory of the dynamics of the 
LV tensor fields themselves, at least at the level of low-energy effective field theory. The most straightforward 
way to construct such a theory is to follow the successful paradigm of scalar field theories with spontaneous 
symmetry breaking, by introducing a tensor potential that is minimized at some non-zero expectation value, 
in addition to a kinetic term for the fields. (Alternatively, it can be a derivative of the field that obtains an 
expectation value, as in ghost condensation models [T5J [121 12D]-) As an additional simplification, we may 
consider models in which the nonzero expectation value is enforced by a Lagrange multiplier constraint, 
rather than by dynamically minimizing a potential; this removes the "longitudinal" mode of the tensor from 
consideration, and may be thought of as a limit of the potential as the mass near the minimum is taken to 
infinity. In that case, there will be a vacuum manifold of zero-energy tensor configurations, specified by the 
constraint. 

All such models must confront the tricky question of stability. Ultimately, stability problems stem from the 
basic fact that the metric has an indefinite signature in a Lorentzian spacetime. Unlike in the case of scalar 
fields, for tensors it is necessary to use the spacetime metric to define both the kinetic and potential terms 
for the fields. A generic choice of potential would have field directions in which the energy is unbounded 
from below, leading to tachyons, while a generic choice of kinetic term would have modes with negative 
kinetic energies, leading to ghosts. Both phenomena represent instabilities; if the theory has tachyons, small 
perturbations grow exponentially in time at the linearized level, while if the theory has ghosts, nonlinear 
interactions create an unlimited number of positive- and negative-energy excitations [21j . There is no simple 
argument that these unwanted features are necessarily present in any model of LV tensor fields, but the 
question clearly warrants careful study. 

In this paper we revisit the question of the stability of theories of dynamical Lorentz violation, and argue 
that most such theories are unstable. In particular, we examine in detail the case of a vector field with a 
nonvanishing expectation value, known as the "aether" model or a "bumblebee" model. For generic choices 
of kinetic term, it is straightforward to show that the Hamiltonian of such a model is unbounded from below, 
and there exist solutions with bounded initial data that grow exponentially in time. 
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There are three specific choices of kinetic term for which the analysis is more subtle. These are the 
sigma-model kinetic term, 

liK^-\d^,A,d^A\ (1) 

which amounts to a set of four scalar fields defined on a target space with a Minkowski metric; the Maxwell 
kinetic term, 

liK = -\F^.F^'' . (2) 
where = ^^J_Al^ — d^A^^ is familiar from electromagnetism; and what we call the "scalar" kinetic term, 

CK^lid.A^r, (3) 

featuring a single scalar degree of freedom. Our findings may be summarized as follows: 

• The sigma-model Lagrangian with the vector field constrained by a Lagrange multiplier to take on 
a timelike expectation value is the only aether theory for which the Hamiltonian is bounded from 
below in every frame, ensuring stability. In a companion paper, we examine the cosmological behavior 
and observational constraints on this model 22]. If the vector field is spacelike, the Hamiltonian is 
unbounded and the model is unstable. However, if the constraint in the sigma-model theory is replaced 
by a smooth potential, allowing the length-changing mode to become a propagating degree of freedom, 
that mode is necessarily ghostlike (negative kinetic energy) and tachyonic (correct sign mass term), 
and the Hamiltonian is unbounded below, even in the timelike case. It is therefore unclear whether 
models of this form can arise in any full theory. 

• In the Maxwell case, the Hamiltonian is unbounded below; however, a perturbative analysis does not 
reveal any explicit instabilities in the form of tachyons or ghosts. The timelike mode of the vector 
acts as a Lagrange multiplier, and there are fewer propagating degrees of freedom at the linear level (a 
"spin-1" mode propagates, but not a "spin-0" mode). Nevertheless, singularities can arise in evolution 
from generic initial data: for a spacelike vector, for example, the field evolves to a configuration in 
which the fixed-norm constraint cannot be satisfied (or perhaps just to a point where the effective 
field theory breaks down). In the timelike case, a certain subset of initial data is well-behaved, but, 
provided the vector field couples only to conserved currents, the theory reduces precisely to conventional 
electromagnetism, with no observable violations of Lorentz invariance. It is unclear whether there exists 
a subset of initial data that leads to observable violations of Lorentz invariance while avoiding problems 
in smooth time evolution. 

• The scalar case is superficially similar to the Maxwell case, in that the Hamiltonian is unbounded 
below, but a perturbative analysis does not reveal any instabilities. Again, there are fewer degrees of 
freedom at the linear level; in this case, the spin-1 mode does not propagate. There is a scalar degree 
of freedom, but it does not correspond to a propagating mode at the level of perturbation theory (the 
dispersion relation is conventional, but the energy vanishes to quadratic order in the perturbations). 
For the timelike aether field, obstacles arise in the time evolution that are similar to those of a spacelike 
vector in the Maxwell case; for a spacelike aether field with a scalar action, the behavior is less clear. 

• For any other choice of kinetic term, aether theories are always unstable. 

Interestingly, these three choices of aether dynamics are precisely those for which there is a unique propagation 
speed for all dynamical modes; this is the same condition required to ensure that the Generalized Second 
Law is respected by a Lorentz- violating theory [231 HI] • 

One reason why our findings concerning stability seem more restrictive than those of some previous analyses 
is that we insist on perturbative stability in all Lorentz frames, which is necessary in theories where the form 
of the Hamiltonian is frame-dependent. In a Lorentz-invariant field theory, it suffices to pick a Lorentz frame 
and examine the behavior of small fluctuations; if they grow exponentially, the model is unstable, while if 
they oscillate, the model is stable. In Lorentz-violating theories, in contrast, such an analysis might miss 
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an instability in one frame that is manifest at the hnear level in some other frame |15l 1251 126| . This can be 
traced to the fact that a perturbation that is "small" in one frame (the value of the perturbation is bounded 
everywhere along some initial spacelike slice), but grows exponentially with time as measured in that frame, 
will appear "large" (unbounded on every spacelike slice) in some other frame. 

As an explicit example, consider a model of a time like vector with a background configuration 

(m, 0,0,0), and perturbations Sa^ = gMe~*"*e'*'''^, where e'' is some constant polarization vector. In this 
frame, we will see that the dispersion relation takes the form 

a;2 = «2p . (4) 

Clearly, the frequency lu will be real for every real wave vector k, and such modes simply oscillate rather 
than growing in time. It is tempting to conclude that models of this form are perturbatively stable for any 
value of V. However, we will see below that when v > I, there exist other frames (boosted with respect to 
the original) in which k can be real but lu is necessarily complex, indicating an instability. These correspond 
to wave vectors for which, evaluated in the original frame, both lu and k are complex. Modes with complex 
spatial wave vectors are not considered to be "perturbations," since the fields blow up at spatial infinity. 
However, in the presence of Lorentz violation, a complex spatial wave vector in one frame may correspond 
to a real spatial wave vector in a boosted frame. We will show that instabilities can arise from initial data 
defined on a constant-time hypersurface (in a boosted frame) constructed solely from modes with real spatial 
wave vectors. Such modes are bounded at spatial infinity (in that frame), and could be superimposed to 
form wave packets with compact support. Since the notion of stability is not frame dependent, the existence 
of at least one such frame indicates that the theory is unstable, even if there is no linear instability in the 
aether rest frame. 

Several prior investigations have considered the question of stability in theories with LV vector fields. Lim 
[9] calculated the Hamiltonian for small perturbations around a constant timelike vector field in the rest 
frame, and derived restrictions on the coefficients of the kinetic terms. Bluhm et al. ^27j also examined 
the timelike case with a Lagrange multiplier constraint, and showed that the Maxwell kinetic term led to 
stable dynamics on a certain branch of the solution space if the vector was coupled to a conserved current. 
It was also found, in [27j, that most LV vector field theories have Hamiltonians that are unbounded below. 
Boundedness of the Hamiltonian was also considered in [21] . In the context of effective field theory, Gripaios 
[55] analyzed small fluctuations of LV vector fields about a flat background. Dulaney, Gresham and Wise [ TT] 
showed that only the Maxwell choice was stable to small perturbations in the spacelike case assuming the 
energy of the linearized modes was non-zero.^ Elliot, Moore, and Stoica [H] showed that the sigma-model 
kinetic term is stable in the presence of a constraint, but not with a potential. 

In the next section, we define notation and fully specify the models we are considering. We then turn to 
an analysis of the Hamiltonians for such models, and show that they are always unbounded below unless the 
kinetic term takes on the sigma-model form and the vector field is timelike. This result does not by itself 
indicate an instability, as there may not be any dynamical degree of freedom that actually evolves along the 
unstable direction. Therefore, in the following section we look carefully at linear stability around constant 
configurations, and isolate modes that grow exponentially with time. In the section after that we show that 
the models that are not already unstable at the linear level end up having ghosts, with the exception of the 
Maxwell and scalar cases. We then examine some features of those two theories in particular. 



II. MODELS 



We will consider a dynamical vector field propagating in Minkowski spacetime with signature ( — I — I — [-). 
The action takes the form 



Sa = J d^x {Ck + Cv) , (5) 



^ This effectively eliminates the scalar case. 
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where Ck is the kinetic Lagrange density and Cv is (minus) the potential. A general kinetic term that is 
quadratic in derivatives of the field can be written^ 

In flat spacetime, setting the fields to constant values at infinity, we can integrate by parts to write an 
equivalent Lagrange density as 



1 A^^ A'^ 

Ck = --/^iF^.i^^'' - P^id^A^^f - Pi^^{d^Ap){d,AP) , (8) 



where F^^^, — df^Ai, — d^A^^ and we have defined 

/3* = /3i + /32 + /33 • (9) 

In terms of these variables, the models specified above with no linear instabilities or negative-energy ghosts 
axe: 

• Sigma model: /3i = /?* , 

• Maxwell: = 0, and 

• Scalar: (3i = 0, 

in all cases with — 0. 

The vector field will obtain a nonvanishing vacuum expectation value from the potential. For most of the 
paper we will take the potential to be a Lagrange multipler constraint that strictly fixes the norm of the 
vector: 

Cv = A(A^A^ ± m^) , (10) 

where A is a Lagrange multiplier whose variation enforces the constraint 

A^A^ = Tm^ . (11) 

If the upper sign is chosen, the vector will be timelike, and it will be spacelike for the lower sign. Later 
we will examine how things change when the constraint is replaced by a smooth potential of the form 
Cv — ^V{Af^) oc ^{Af^A^ ± m^)^. It will turn out that the theory defined with a smooth potential is only 
stable in the limit as ^ — > c». In any case, unless we specify otherwise, we assume that the norm of the 
vector is determined by the constraint (111. 
We are left with an action 



Sa = d X 



1 A>^ A 

-^PiF^.F>^- - P.id^A^'f - /34 :^{d^Ap){d,AP) + \{A>^A^ ± m") 



(12) 



The Euler-Lagrange equation obtained by varying with respect to A^ is 

Pl^^F^''' + M^'d^A^' + piG" - -A^" , (13) 

where we have defined 

= ^ [A^{d>,A'')F,- + A''{dxA^d,A'' + A^d^^d^A^)] . (14) 
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In terms of the coefficients, Ci, defined in [8] and used in many other publications on asther theories, 



where G is the gravitational constant. 
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Since the fixed-norm condition (111 is a constraint, we can consistently plug it back into the equations of 
motion. Multiplying (13 1 by A^, and using the constraint, we can solve for the Lagrange multiplier. 
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(15) 



Inserting this back into ( 13 ), we can write the equation of motion as a system of three independent equations: 

ApAij 



Qp = [iipu± 



(16) 



The tensor r]p,jzLm~'^ApAi, acts to take what would be the equation of motion in the absence of the constraint, 
and project it into the hyperplane orthogonal to A^^. There are only three independent equations because 
APQp vanishes identically, given the fixed norm constraint. 



A. Validity of effective field theory 

As in this paper we will restrict our attention to classical field theory, it is important to check that any 
purported instabilities axe found in a regime where a low-energy effective field theory should be valid. The 
low-energy degrees of freedom in our models are Goldstone bosons resulting from the breaking of Lorentz 
invariance. The effective Lagrangian will consist of an infinite series of terms of progressively higher order 
in derivatives of the fields, suppressed by appropriate powers of some ultraviolet mass scale M. If we were 
dealing with the theory of a scalar field <&, the low-energy effective theory would be valid when the canonical 
kinetic term (9$)^ was large compared to a higher-derivative term such as 

For fluctuations with wavevector k'^ — (uj,k), we have ~ k^, and the lowest-order terms accurately 
describe the dynamics whenever |fc| < M. A fluctuation that has a low momentum in one frame can, of 
course, have a high momentum in some other frame, but the converse is also true; the set of perturbations 
that can be safely considered "low-energy" looks the same in any frame. 

With a Lorentz- violating vector field, the situation is altered. In addition to higher-derivative terms of the 
form M~^(9^A)^, the possibility of extra factors of the vector expectation value leads us to consider terms 
such as 

C^^A'id'Af. (18) 

The number of such higher dimension operators in the effective field theory is greatly reduced because 
Ap^A'^ = — and, therefore, Apd^.A'^ = 0. It can be shown that an independent operator with n derivatives 
includes at most 2n vector fields, so that the term highlighted here has the largest number of A's with four 
derivatives. We expect that the ultraviolet cutoff M is of order the vector norm, M k, m. Hence, when we 
consider a background timelike vector field in its rest frame, 

= (m, 0,0,0), (19) 

the £4 term reduces to m^'^ {d'^ A)'^ , and the effective field theory is valid for modes with k < m, just as in 
the scalar case. 

But now consider a highly boosted frame, with 

= ('T^cosh?7, msinh?], 0, 0) . (20) 

At large 77, individual components of A will scale as el''' , and the higher-derivative term schematically becomes 

£4 ~ ^eS'"'(9'A)2. (21) 
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For modes with spatial wave vector fc = |fc| (as measured in this boosted frame), we are therefore comparing 
^-2g6|r(|^4 -^j^jj canonical term k^. The lowest-order terms therefore only dominate for wave vectors 
with 

fc<e-3|''lm. (22) 

In the presence of Lorentz violation, therefore, the realm of validity of the effective field theory may be 
considerably diminished in highly boosted frames. We will be careful in what follows to restrict our con- 
clusions to those that can be reached by only considering perturbations that are accurately described by 
the two-derivative terms. The instabilities we uncover are infrared phenomena, which cannot be cured by 
changing the behavior of the theory in the ultraviolet. We have been careful to include all of the lowest order 
terms in the effective field theory expansion — the terms in ([8]). 



III. BOUNDEDNESS OF THE HAMILTONIAN 



We would like to establish whether there are any values of the parameters /3i , /3* and /34 for which the 
aether model described above is physically reasonable. In practice, we take this to mean that there exist 
background configurations that are stable under small perturbations. It seems hard to justify taking an 
unstable background as a starting point for phenomenological investigations of experimental constraints, as 
we would expect the field to evolve on microscopic timescales away from its starting point. 

"Stability" of a background solution Xq to a set of classical equations of motion means that, for any 
small neighborhood Uq of Xq in the phase space, there is another neighborhood Ui of Xq such that the time 
evolution of any point in Uq remains in Ui for all times. More informally, small perturbations oscillate around 
the original background, rather than growing with time. A standard way of demonstrating stability is to 
show that the Hamiltonian is a local minimum at the background under consideration. Since the Hamiltonian 
is conserved under time evolution, the allowed evolution of a small perturbation will be bounded to a small 
neighborhood of that minimum, ensuring stability. Note that the converse does not necessarily hold; the 
presence of other conserved quantities can be enough to ensure stability even if the Hamiltonian is not 
bounded from below. 

One might worry about invoking the Hamiltonian in a theory where Lorentz invariance has been sponta- 
neously violated. Indeed, as we shall see, the form of the Hamiltonian for small perturbations will depend 
on the Lorentz frame in which they are expressed. To search for possible linear instabilities, it is necessary 
to consider the behavior of small perturbations in every Lorentz frame. 



The Hamiltonian density, derived from the action (12) via a Legendre transformation, is 



A3 Ak /lO/lO 

+ PA^{djAp){dkAP) - /34^(aoAp)(5oA^), (24) 
where Latin indices i,i run over {1,2,3}. The total Hamiltonian corresponding to this density is 
H = ( (fixTL 

(fx{l3i{d^A,d^,A, - a^^o^^^o) + (/3i - /3*)[(aoAo)2 - [d.A^f] 

+ P^^^{d,Ap){dkAP) ~ p^^^{doAp){d^AP)) . (25) 

We have integrated by parts and assumed that diAj vanishes at spatial infinity; repeated lowered indices 
are summed (without any factors of the metric). Note that this Hamiltonian is identical to that of a theory 
with a smooth (positive semi-definite) potential instead of a Lagrange multiplier term, evaluated at field 
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configurations for which the potential is minimized. Therefore, if the Hamiltonian is unbounded when the 
fixed-norm constraint is enforced by a Lagrange multipHer, it will also be unbounded in the case of a smooth 
potential. 

There are only three dynamical degrees of freedom, so we may reparameterize such that the fixed- norm 
constraint is automatically enforced and the allowed three-dimensional subspace is manifest. We define a 
boost variable (j) ^^nd angular variables 6 and tp, so that we can write 

= TO cosh (f) (26) 
A, = TOsinh0/,(0,V') (27) 



in the timelike case with A^A^ ~ 



AQ = msm\\(j) (28) 
Ai = m cosh cj)fi{e,^) (29) 



in the spacelike case with A^A^^ — +m? . In these expressions, 

/i = cos 6* cos i/' (30) 
/2 = cosOsmip (31) 
h = sine, (32) 

so that fifi — 1. In terms of this parameterization, the Hamiltonian density for a timelike aether field 
becomes 

— ^ = /3i sinh^ cj)dj,d^f, + /3i 9^09^0 + (/3i - /?*) [(^o^)' sinh^ - (cosh + sinh ^dj,f] 

+ /34 sinh^ cf, [(/,9.0)2 + sinh^ cl>{fAfi){fjd,fi)] - (3^ cosh^ [(^o^)" + sinh^ 0(9o/»)'] , (33) 
while for the spacelike case we have 

= Pi cosh^ (l^d^Jid^fi - Pidf,(t)df,(t) + (/3i - /?*) [(9o0)^ cosh^ (/j - (sinh (j)fid^(t) + cosh (t)difi)^] 

- Acosh^ [(/.a,0)2 - cosh^ cl>{fAfi){fjdjfi)] + /34sinh2 [(900)^ - cosh^ 0(9o/^)'] ■ (34) 

Expressed in terms of the variables 0, 9, ip, the Hamiltonian is a function of initial data that automatically 
respects the fixed-norm constraint. We assume that the derivatives dfj_A^(tQ, x) vanish at spatial infinity. 



A. Timelike vector field 



We can now determine which values of the parameters /34} lead to Hamiltonians that are bounded 

below, starting with the case of a timelike aether field. We can examine the various possible cases in turn. 

• Case One: f3i — /3* and /S^ — 0. 

This is the sigma- model kinetic term ([T]). In this case the Hamiltonian density simplifies to 

= TO2/3i(sinh2 cl)d^f,dj, + d^<i>d^(j>) . (35) 

It is manifestly non-negative when (3i > 0, and non-positive when /3i < 0. The sigma-model choice 
/3i = /9* > therefore results in a theory that is stable. (See also §6.2 of [10 .) 

• Case Two: /3i < and /34 = 0. 

In this case, consider configurations with (9o/i) 7^ 0, (difj) = 0, d^cj) — 0, sinh^ 0^1. Then we have 

- m^Pi sinh^ <t>{daf,f. (36) 
For /3i < 0, the Hamiltonian can be arbitrarily negative for any value of 
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FIG. 1: Hamiltonian density (vertical axis) when l3i — 1, p, = 1.1, and 9 = il) = dycf) = dz4> = as a function of dt4> 
(axis pointing into page) and d^cf) (axis pointing out of page) for various (j> ranging from zero to (jicrit = tanh~^ ^ /3i//3t, 
the value of 4> for which the Hamiltonian is flat at dx4> = 0, and beyond. Notice that the Hamiltonian density turns 
over and becomes negative in the dt4' direction when (j) > ^crit. 



• Case Three: /3i > 0, < /3i, and {3^ = 0. 

We consider configurations with d^fi = 0, fidi(t) 7^ 0, do(f> = 0, cosh^ </) ^ 1, which gives 

- m2(/3, - /3i) cosh^ ^ifA<pf. (37) 
Again, this can be arbitrarily negative. 

• Case Four: Pi > 0, /5* > /3i, and Pi = 0. 

Now we consider configurations with df^ft = 0, fidicf) = 0, do4> ^ 0, sinh^ ^ ^ 1. Then, 

~ m^iPi - /3*) sinh^ (/)(5o0)^ (38) 

which can be arbitrarily negative. 

• Case Five: ^4 7^ 0. 

Now we consider configurations with 7^ 0, = and sinh^ (j) ^ I. Then, 

H(') ~ m2/?4 [sinh^ <l>{f.Afi){fkdkfi) - sinh^ ^cosh^ </>(ao/,)'] , (39) 

which can be arbitrarily negative for any non-zero (3^ and for any values of /3i and 

For any case other than the sigma- model choice /3i = /?* , it is therefore straightforward to find configurations 
with arbitrarily negative values of the Hamiltonian. 

Nevertheless, a perturbative analysis of the Hamiltonian would not necessarily discover that it was un- 
bounded. The reason for this is shown in Fig. [T| which shows the Hamiltonian density for the theory with 
Pi = 1, /3* = 1.1, in a restricted subspace where dycf) = dz4> = and 6 = cj) = 0, leaving only (j>, dtcj), and 
dx4> as independent variables. We have plotted 7i as a function of dt<j> and dx<p for four different values of (f>. 
When (j) is sufficiently small, so that the vector is close to being purely timelike, the point dttp — dx4> = is a 
local minimum. Consequently, perturbations about constant configurations with small (j) would appear sta- 
ble. But for large values of 0, the unboundedness of the Hamiltonian becomes apparent. This phenomenon 
will arise again when we consider the evolution of small perturbations in the next section. At the end of this 
section, we will explain why such regions of large are still in the regime of validity of the effective field 
theory expansion. 
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B. Spacelike vector field 

We now perform an equivalent analysis for an aether field with a spacelike expectation value. In this case 



all of the possibiUties lead to Hamiltonians (34 1 that are unbounded below, and the case /3i = /?, > is not 
picked out. 

• Case One: /3i < and /34 = 0. 

Taking (9^(/)) = 0, d,f, = 0, dol, ^ 0, we find 

H*"' - m'^Pi cosh^ <p{dof,f. (40) 

• Case Two: /3i > 0, /3* < /3i, and Pi = 0. 

Now we consider d^fi — Q, di(f> ^ 0, do(f> — 0, giving 

^ ^2 + (^_^ _ 0(/,5,0)2] . (41) 

• Case Three: /3i > 0, /3* > and /34 = 0. 

In this case we examine (docj)) ^ 0, 9^/^ = 0, di(\> — 0, which leads to 

n^''^ - m2(/5i - /3,) cosh^ (/'(ao0)^- (42) 

• Case Four: fi^ ^ 0. 

Now we consider configurations with df^fi^Q,df^(j) = Q and sinh^ (f> ^ 1. Then, 

^ m2/34 (cosh^ Hmmfkdkfi) - cosh^ ^sinh^ 0(ao/^)') ■ (43) 

In every case, it is clear that we can find initial data for a spacelike vector field that makes the Hamiltonian 
as negative as we please, for all possible P4 and 

C. Smooth Potential 

The usual interpretation of a Lagrange multiplier constraint is that it is the low-energy limit of smooth 
potentials when the massive degrees of freedom associated with excitations away from the minimum cannot 
be excited. We now investigate whether these degrees of freedom can destabilize the theory. Consider the 
most general, dimension four, positive semi-definite smooth potential that has a minimum when the vector 
field takes a timelike vacuum expectation value. 



V^^-iA,A>^ + m^)^, (44) 



where ^ is a positive dimensionless parameter. The precise form of the potential should not affect the results 
as long as the potential is non-negative and has the global minimum at A/^^A^ = —m^. 

We have seen that the Hamiltonian is unbounded from below unless the kinetic term takes the sigma-model 
form, {df^A^){d^A'^). Thus we take the Lagrangian to be 



£ = -\{d,A,){d^A'^) - |(A,A^ + m'f. (45) 



Consider some fixed timelike vacuum A^ satisfying Af^A'^ = —m?. We may decompose the aether field 
into a scaling of the norm, represented by a scalar and an orthogonal displacement, represented by vector 
satisfying A^B'^ — 0. We thus have 

A^ = A^~^ + B,,, (46) 
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where 



B^=[v^,+ ^)a'' and ^^^^+m. (47) 



With this parameterization, the Lagrangian is 

£ = J(a^<i>)(9^<i>) - ^{dMid^B'') - |(2m<i> + S^B^ - ^^)^. (48) 



The field $ automatically has a wrong sign kinetic term, and, at the linear level, propagates with a dispersion 
relation of the form 

ujI = P - 2^m^. (49) 

We see that in the case of a smooth potential, there exists a ghostlike mode (wrong-sign kinetic term) that is 
also tachyonic with spacelike wave vector and a group velocity that generically exceeds the speed of light. It 
is easy to see that sufficiently long-wavelength perturbations will exhibit exponential growth. The existence 
of a ghost when the norm of the vector field is not strictly fixed was shown in fl4J . 

In the limit as ^ goes to infinity, the equations of motion enforce a fixed-norm constraint and the ghostlike 
and tachyonic degree of freedom freezes. The theory is equivalent to one of a Lagrange multiplier if the limit 
is taken appropriately. 



D. Discussion 



To summarize, we have found that the action in ( [T2| leads to a Hamiltonian that is globally bounded from 
below only in the case of a timelike sigma-model Lagrangian, corresponding to /3i = /?* > and /?4 = 0. 
Furthermore, we have verified (as was shown in [H]) that if the Lagrange multiplier term is replaced by a 
smooth, positive semi-definite potential, then a tachyonic ghost propagates and the theory is destabilized. 

If the Hamiltonian is bounded below, the theory is stable, but the converse is not necessarily true. The 
sigma-model theory is the only one for which this criterion suffices to guarantee stability. In the next section, 
we will examine the linear stability of these models by considering the growth of perturbations. Although 
some models are stable at the linear level, we will see in the following section that most of these have 
negative-energy ghosts, and are therefore unstable once interactions are included. The only exceptions, both 
ghost-free and linearly stable, are the Maxwell Q and scalar ^ models. 

We showed in the previous section that, unless /3* — f3i and (34 are exactly zero, the Hamiltonian is 
unbounded from below. However, the effective field theory breaks down before arbitrarily negative values 
of the Hamiltonian can be reached; when (3^ ^ Pi and/or P4 ^ 0, in regions of phase space in which Ti < 
(schematically) , 

-TO2e-*l'^l(ae)2 where Q e {<j>,9,il;}. (50) 

The effective field theory breaks down when kinetic terms with four derivatives (the terms of next highest 
order in the effective field theory expansion) are on the order of terms with two derivatives, or, in the angle 
parameterization, when 

TO2e4|<^l(ae)2~e«l*l(9e)^ (51) 

In other words, the effective field theory is only valid when 

e2|"*l|ae| < m. (52) 

In principle, terms in the effective action with four or more derivatives could add positive contributions to 
the Hamiltonian to make it bounded from below. However, our analysis shows that the Hamiltonian (in 
models other than the timelike sigma model with fixed norm) is necessarily concave down around the set of 
configurations with constant aether fields. If higher-derivative terms intervene to stabilize the Hamiltonian, 
the true vacuum would not have iJ = 0. Theories could also be deemed stable if there are additional 
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symmetries that lead to conserved currents (other than energy-momentum density) or to a reduced number 
of physical degrees of freedom. 

Regardless of the presence of terms beyond leading order in the effective field theory expansion, due to the 
presence of the ghost-like and tachyonic mode (found in the previous section), there is an unavoidable problem 
with perturbations when the field moves in a smooth, positive semi-definite potential. This exponential 
instability will be present regardless of higher order terms in the effective field theory expansion because it 
occurs for very long- wavelength modes (at least around constant-field backgrounds). 



IV. LINEAR INSTABILITIES 



We have found that the Hamiltonian of a generic aether model is unbounded below. In this section, we 
investigate whether there exist actual physical instabilities at the linear level — i.e., whether small pertur- 
bations grow exponentially with time. It will be necessary to consider the behavior of small fluctuations in 
every Lorentz frame, not only in the aether rest frame |15l 125} I26j . We find a range of parameters Pi for 
which the theories are tachyon-free; these correspond (unsurprisingly) to dispersion relations for which the 
phase velocity satisfies < < 1. In fjvlwe consider the existence of ghosts. 



A. Timelike vector field 



Suppose Lorentz invariance is spontaneously broken so that there is a preferred rest frame, and imagine 
that perturbations of some field in that frame have the following dispersion relation: 



V '^uP' — k ■ k. 



This can be written in frame-invariant notation as 



(53) 



(54) 



where is a timelike Lorentz vector that characterizes the 4-velocity of the preferred rest frame. So, in the 
rest frame, — {1,0,0,0}. Indeed, in the Appendix, we find dispersion relations for the aether modes of 
exactly the form in (54 1 with — /m and (A27l 



(3i 



Pi -Pa 



and (A28I 



P* 



Pi- Pi' 



(55) 



(56) 



Now consider the dispersion relation for perturbations of the field in another ("primed") frame. Let's solve 
for fcg = Lo' , the frequency of perturbations in the new frame. Expanded out, the dispersion relation reads 



Lo'^{l + {v-^ - 1)(<'°)2) + 2w'(«-2 _ i)t'h"k[ -k'-k' + {v~^ - l){t''k[f = 
where i S {1, 2, 3}. The solution for us' is: 



l-H(w-2-l)(t'0)2 



(57) 



(58) 



The theory of perturbations about a constant background is equivalent to a theory with expHcit Lorentz violation because 
the first order Lagrange density includes the term, XA^SAf^, where is effectively some constant coefficient. 
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where 



D 



= k'.k'+ {v-^ - 1) . k' - (f"fc^)2) 



(59) 



In general, t'" — coshrj and t'* = sinhr^n', where fiin' = 1 and 77 = cosh ^ 7 is a boost parameter. We 
therefore have 



(60) 



where fc' = k'/\k'\. Thus D(t) is clearly greater than zero if u < 1. However, if u > 1 then Z3(f) can be 

negative for very large boosts if k' is not parallel to the boost direction. 

The sign of the discriminant determines whether the frequency lo' is real- or complex- valued. We 
have shown that when the phase velocity v of some field excitation is greater than the speed of light in a 
preferred rest frame, then there is a (highly boosted) frame in which the excitation looks unstable — that is, 
the frequency of the field excitation can be imaginary. More specifically, plane waves traveling along the 
boost direction with boost parameter 7 = coshry have a growing amplitude if 7^ > 1/(1 — u~^) > 0. 

In Appendix |Aj we find dispersion relations of the form in (54 1 for the various massless excitations about a 
constant timelike background {t^ — A'^/m). Requiring stability and thus < f ^ < 1 leads to the inequalities, 



< 



Pi 
Pi -Pi 



< 1 



and 



0< ^* < 1 

- Pi -Pi- 



rn 



(62) 



Models satisfying these relations are stable with respect to linear perturbations in any Lorentz frame. 



B. Spacelike vector field 



We show in Appendix [X] that fluctuations about a spacelike, fixed-norm, vector field background have 
dispersion relations of the form 



with = A^'/m and (A27l 



and (A28I 



{v^ - l)(s'^fc^)2 = -fc^fc^ 

2 Pi + Pi 
^ = — 

2 P1 + P4 



(63) 



(64) 



(65) 



In frames where — {0, s}, v is the phase velocity in the s direction. 



Consider solving for fcp = w' in an arbitrary ("primed") frame. The solution is as in (58 1, but with 



and t''^ -> s'^. Thus, 



OJ 



where 



l + (l-'y2)(s'0)2 
D(,) ^k' -k'^ {v^ - 1) \{s"yk' ■ k' - (s"fcO' 



(66) 



(67) 



13 



In general, s"^ — sinhyy and s" = cosh 77 where n^n' = 1 and 77 — cosh ^ 7 is a boost parameter. So, 



which can be rewritten, 



D(^) = fc' • fc' |l - (u^ - 1) sinh^ 7? - cosh^ 77 (n • k'f | 
= fc' • fc' ji-^ + (1 - v'^) cosh^ 77 1 - (n • fc')^ } • 



(68) 



(69) 



It is clear that -D(s) is non-negative for all values of 77 if and only if < < 1. The theory will be unstable 
unless < < 1. 

The dispersion relations of the form ( 63 1 for the massless excitations about the spacelike background are 
given in Appendix |A] The requirement that < < 1 implies 



0<^^<1 



(70) 



and 



0<^<1, 



(71) 



Models of spacelike aether fields will only be stable with respect to linear perturbations if these relations are 
statisfied. 

The requirements (62) or (711 do not apply in the Maxwell case (when /3* = = /S^), and those of (61 1 or 



(70 1 do not apply in the scalar case (when /3i = = Pi), since the corresponding degrees of freedom in each 
case do not propagate. 



C. Stability is not frame-dependent 

The excitations about a constant background are massless {i.e. the frequency is proportional to the mag- 
nitude of the spatial wave vector), but they generally do not propagate along the light cone. In fact, when 
V > 1, the wave vector is timelike even though the cone along which excitations propagate is strictly outside 
the light cone. We have shown that such excitations blow up in some frame. The exponential instability 
occurs for observers in boosted frames. In these frames, portions of constant-time hypersurfaces are actually 
inside the cone along which excitations propagate. 

Why do we see the instability in only some frames when performing a linear stability analysis? Consider 
boosting the wave four-vectors of such excitations with complex-valued frequencies and real-valued spatial 
wave vectors back to the rest frame. Then, in the rest frame, both the frequency and the spatial wave 
vector will have non-zero imaginary parts. Such solutions with complex-valued k require initial data that 
grow at spatial infinity and are therefore not really "perturbations" of the background. But even though the 
asther field defines a rest frame, there is no restriction against considering small perturbations defined on a 
constant-time hypersurface in any frame. Well-behaved initial data can be decomposed into modes with real 
spatial wave vectors; if any such modes lead to runaway growth, the theory is unstable. 



V. NEGATIVE ENERGY MODES 



We found above that manifest perturbative stability in all frames requires < < 1- In the Appendix, 
we show that there are two kinds of propagating modes, except when (3^, — (34 — or when (3i = P4 — 0. 
Based on the dispersion relations for these modes, the < < 1 stability requirements translated into the 
inequalities for /3*,/3i, and P4 in (61 1- (62) for timelike aether and (70)- (71 1 for spacelike aether. We shall 
henceforth assume that these inequalities hold and, therefore, that ui and k for each mode are real in every 
frame. We will now show that, even when these requirements are satisfied and the theories are linearly 
stable, there will be negative-energy ghosts that imply instabilities at the nonlinear level (except for the 
sigma model. Maxwell, and scalar cases). 
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For timelike vector fields, with respect to the aether rest frame, the various modes correspond to two spin-1 
degrees of freedom and one spin-0 degree of freedom. Based on their similarity in form to the timelike aether 
rest frame modes, we will label these modes once and for all as "spin-1" or "spin-0," even though these 
classifications are only technically correct for timelike fields in the aether rest frame. 

The solutions to the first order equations of motion for perturbations 6A^ about an arbitrary, constant, 
background satisfying A^Afj_±m'^ =0 are (see Appendix [a| : 



d'kq^{k)e^'^^^\ q,{k)^q*J-k) 



where either, 



q^{k)^ia''kP — e^,p, and p^k^k'' + 13^-^^^^ 



TO 



and a" Ay = (spin-1) 



(72) 



(73) 



where a'^ are real-valued constants or 
A,.A 



k" 



and 



/3*?7m- + (/34±(/3*-/3i)) 



k^'k" = (spin-0) (74) 



where a is a real- valued constant. 

Note that when fii — = Q, corresponding to the scalar form of (|3|, the spin-1 dispersion relation is 
satisfied trivially, because the spin-1 mode does not propagate in this case. Similarly, when /3* = /34 = 0, 
the kinetic term takes on the Maxwell form in ([2]) and the spin-0 dispersion relation becomes A^k'^ — 0; the 
spin-0 mode does not propagate in that case. 



The Hamiltonian (p5|) for cither of these modes is 
(fk 



H 



{ [/3i(l^2 + ^ . fc) + p^{-(fu;f + {a'k,f)\ q^^q^ + (/?i - /?*) (c^^g^^O + hqlk^q,)] 



(75) 



where k^ — lo — uj{k) is given by the solution to a dispersion relation and where = A^ jm. One can show 
that, as long as /3i and (5/^ satisfy the conditions (61 1 or (70 1 that guarantee real frequencies w in all frames, 
we will have 



(76) 



for all timelike and spacelike vector perturbations. We will now proceed to evaluate the Hamiltonian for 
each mode in different theories. 



A. Spin-1 energies 



In this section we consider nonvanishing /34, and show that the spin-1 mode can carry negative energy 
even when the conditions for linear stability are satisfied. 
Timelike vector field. Without loss of generality, set 



A^ = TO (cosh ?7, sinh 77 h). 
where fi ■ n = 1. The energy of the spin-1 mode in the timelike case is given by 

2X =F /34sinh(277)(n • k)^/X 



H = / d'k{k ■ k)q:q^ 



Pi — Pi cosh J] 



where 



X 



/3i |/3i + Pa: {fi-kf sinh^ - cosh^ 77 | 



(77) 



(78) 



(79) 
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Looking specifically at modes for which n • fc = +1, we find 

" 2/?i(/3i - Pi) T /34sinh(277) - pj) 

Pi — Pi cosh^ r] 



H = / d'k{k ■ k)q*X 



(80) 



The energy of such a spin-1 perturbation can be negative when |/94 sinh(2r7)| > 2^J Pi{Pi — Pi). Thus it 
is possible to have negative energy perturbations whenever Pi ^ Q. Perturbations with wave numbers 
perpendicular to the boost direction have positive semi-definite energies. 

Spacelike vector field. Without loss of generality, for the spacelike case we set 



A., 



m(sinh 77, cosh -q fi) 



where fi ■ n = 1. The energy of the spin-1 mode in this case is given by 



H 



d^'kik ■ k)q;q^ 



2X =F ^4sinh(277)(ri • k)VX 
Pi — Pi sinh^ 77 



where 



X ^ Pi ^Pi + Pi (n • fc)^ cosh^ 77 - sinh^ r; 
Looking at modes for which n • = -1-1, we find 

d^k{k ■ k)qlq^' 



H 



2Pi{Pi + Pi) T /?4sinh(277)v//3i(/3i + /34) 



Pi — Pi sinh ?7 



(81) 



(82) 



(83) 



(84) 



Thus, the energy of perturbations can be negative when |/34 sinh(2?7)| > 2yJ Pi{Pi + Pi). Thus it is possible 
to have negative energy perturbations whenever Pi^ Q. Perturbations with wave numbers perpendicular to 
the boost direction have positive semi-definite energies. In either the timelike or spacelike case, models with 
Pi^Q feature spin-1 modes that can be ghostlike. 

We note that the effective field theory is valid when k < e^'^l'^'m, as detailed in {II A But even if rj is very 
large, the effective field theory is still valid for very long wavelength perturbations, and therefore such long 
wavelength modes with negative energies lead to genuine instabilities. 



B. Spin-0 energies 



We now assume the inequalities required for linear stability, (62 1 or (71 1, and also that Pi = 0. We showed 
above that, otherwise, there are growing modes in some frame or there are propagating spin-1 modes that 
have negative energy in some frame. When /3* 7^ 0, the energy of the spin-0 mode in (74 1 is given by 

H = 2Pia^ J d^k (apkP)^ (u!^{k) [±1 - (1 - Pi/P*)al] + Lo{k) ao(l - PilP*)a,k,^ (85) 

for Afj_A^ ± — and = Ap_/m. 

Timelike vector field. We will now show that the quadratic order Hamiltonian can be negative when the 
background is timelike and the kinetic term does not take one of the special forms (sigma model. Maxwell, 
or scalar). Without loss of generality we set ao — cosh 77 and cii = sinh77?T,i, where n ■ h — 1. Then plugging 
the freqency oj{k), as defined by the spin-0 dispersion relation, into the Hamiltonian (|85[) gives 



H = Pia^ / d^kiapkPy^ 



where 



2X±{1~ Pi/P^) sinh2r;(ri • k)VX 
l + {Pi/P, - l)cosh2r7 

X = l + {Pi IP* - l)[cosh2 77 - (n • fc)2 sinh^ 77]. 



(86) 



(87) 
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If n • /c ^ 0, the energy can be negative. In particular, if n • fc = 1 we have 

"2/3i//3, ± (1 - /3i//3,) sinh277yAM 



l + (/3i//3* - l)cosh"?7 



Given that — 1 > 0, can be negative when | sinh2?7| > 2^ — !)• 

We have thus shown that, for timeHke backgrounds, there are modes that in some frame have negative 

energies and/or growing amphtudes as long as /3i 7^ /3i ^ 0, and /3* 7^ 0. Therefore, the only possibly 

stable theories of timelike aether fields are the special cases mentioned earlier: the sigma-model {(3i = (3^,), 

Maxwell (/3* = 0), and scalar {Pi — 0) kinetic terms. 

Spacelike vector field. For the spacelike case, without loss of generality we set So = sinhyy and — 

coshrihi, where h - h — 1. Once again, plugging the frequency Lu{k) into the Hamiltonian (85) gives 



where 



± (1 - /3i//3*) sinh 277(n • k)y/X 
l + (l-/3i//3*)sinh^7; 



X = l + {1- sinh^ r]-{fi- kf cosh^ r? 



(89) 



(90) 



Upon inspection, one can see that there are values oi h ■ k and rj that make H negative, except when = 
(Maxwell) or Pi — (scalar). Again, the Hamiltonian density is less than zero for modes with wavelengths 
sufficiently long (fc < e"^''''™), so the effective theory is valid. 



VI. MAXWELL AND SCALAR THEORIES 



We have shown that the only version of the aether theory ( 12 1 for which the Hamiltonian is bounded below 
is the timelike sigma-model theory £x = —{^/'^){dfj,A^){d'^A'^, corresponding to the choices Pi = P^ = 0, 
with the fixed-norm condition imposed by a Lagrange multiplier constraint. (Here and below, we rescale the 
field to canonically normalize the kinetic terms.) However, when we looked for explicit instabilities in the 
form of tachyons or ghosts in the last two sections, we found two other models for which such pathologies 
are absent: the Maxwell Lagrangian 



1 



K 



corresponding to = = /?4, and the scalar Lagrangian 



(91) 



(92) 



corresponding to /3i = = /34. In both of these cases, we found that the Hamiltonian is unbounded below,** 
but a configuration with a small positive energy does not appear to run away into an unbounded region of 
phase space characterized by large negative and positive balancing contributions to the total energy. 

These two models are also distinguished in another way: there are fewer than three propagating degrees of 
freedom at first order in perturbations in the Maxwell and scalar Lagrangian cases, while there are three in all 
others. This is closely tied to the absence of perturbative instabilities; the ultimate cause of those instabilities 
can be traced to the difficulty in making all of the degrees of freedom simultaneously well-behaved. The 
drop in number of degrees of freedom stems from the fact that Ao lacks time derivatives in the Maxwell 
Lagrangian and that the Ai lack time derivatives in the scalar Lagrangian. In other words, some of the 
vector components are themselves Lagrange multipliers in these special cases. 



^ Boundedness of the Hamiltonian was considered in [30) . 
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Only two pcrturbative degrees of freedom — the spin-1 modes — propagate in the Maxwell case (cf. (73)- 
(74 1 when = = /?4). The "mode" in (74 1 is a gauge degree of freedom; at first order in perturbations 
the Lagrangian has a gauge-like symmetry under 5Af^ — > 6Af^ + d^(j){x) where A^^d^cj) = 0- expected 
of a gauge degree of freedom, the spin-0 mode has zero energy and does not propagate. Meanwhile, the 
spin-1 perturbations propagate as well-behaved plane waves and have positive energy. We note that the 
Dirac method for counting degrees of freedom in constrained dynamical systems implies that there are three 
degrees of freedom [27].^ The additional degree of freedom, not apparent at the linear level, could conceivably 
cause an instability; this mode does not propagate because it is gauge-like at the linear level, but there is no 
gauge symmetry in the full theory. 

In the scalar case, there are no propagating spin-1 degrees of freedom. The spin-0 degree of freedom has 
a nontrivial dispersion relation but no energy density (cf. (73 1- (74), (86), and (89) when f3i — Q = (5^) at 
leading order in the perturbations. Essentially, the fixed-norm constraint is incompatible with what would 
be a single propagating scalar mode in this model; the theory is still dynamical, but perturbation theory 
fails to capture its dynamical content. 

Each of these models displays some idiosyncratic features, which we now consider in turn. 



A. Maxwell action 



The equation of motion for the Maxwell Lagrangian with a fixed-norm constraint is 

Q^pt,^ = -2\A'' . (93) 
Setting Afj^A^ = ^m?, the Lagrange multiplier is given by 

A = ±2^^.5,F-''- (94) 

For timelike aether fields, the sign of A is preserved along timelike trajectories since, when the kinetic term 
takes the special Maxwell form, there is a conserved current (in addition to energy-momentum density) due 
to the Bianchi identity^: 

o = a,(a^i^^n--2a.(AA^). (95) 

In particular, the condition that A = is conserved along timelike A'^ [5JI17]. In the presence of interactions 
this will continue to be true only if the coupling to external sources takes the form of an interaction with a 
conserved current, A^J^^ with d^J^ = 0. 

If we take the timelike Maxwell theory coupled to a conserved current and restrict to initial data satisfying 
A = at every point in space, the theory reduces precisely to Maxwell electrodynamics — not only in the 
equation of motion, but also in the energy-momentum tensor. We can therefore be confident that this theory, 
restricted to this subset of initial data, is perfectly well-behaved, simply because it is identical to conventional 
electromagnetism in a nonlinear gauge ES| ^37 . 

In the case of a spacelike vector expectation value, there is an explicit obstruction to finding smooth time 
evolution for generic initial data. In this case, the constraint equations are 

~Al + A,A, = and 9,aMo - d^diA' = -2\Aq. (96) 

Suppose spatially homogeneous initial conditions for the A^ are given. Without loss of generality, we can 
align axes such that 

^m(*o) = (^(io), 0,0,^3(^0)), (97) 



^ For a discussion of constrained dynamical systems see 1311 . 

® If A > initially, then it must pass through A = to reach A < — but A = is conserved along timelike trajectories, so A 
can at best stop at A = 0. 
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where — + A| — . If AiAi ^ m? , the equations of motion are 

= 0. (98) 



The J/ = 3 equation reads 



whose solutions are given by 



d,F>^, = -^^Q, (99) 



A:i{t)^ Az{to) + C{t-to), (100) 



where C is determined by initial conditions. Aq is determined by the fixed-norm constraint Aq — 
±^/A^^^rn^ . If C 7^ 0, Aq will eventually evolve to zero. Beyond this point, A3 keeps decreasing, and 
the fixed-norm condition requires that Aq be imaginary, which is unacceptable since A^ is a real-valued 
vector field. Note that this never happens in the timelike case, as there always exists some real Aq that 
satisfies the constraint for any value of A3. The problem is that A^ evolves into the ball Af < m^, which 
is catastrophic for the spacelike, but not the timelike, case. An analogous problem arises even when the 
Lagrange multiplier constraint is replaced by a smooth potential. 

It is possible that this obstruction to a well-defined evolution will be regulated by terms of higher order 
in the effective field theory. Using the fixed-norm constraint and solving for Aq, the derivative is 

d,Ao = ,^ f d,A,. (101) 

As AjAj approaches m^, with finite derivatives of the spatial components, the derivative of the Aq component 
becomes unbounded. If higher-order terms in the effective action have time derivatives of the component 
these terms could become relevant to the vector field's dynamical evolution, indicating that we have left 
the realm of validity of the low-energy effective field theory we are considering. 

We axe left with the question of how to interpret the timelike Maxwell theory with intial data for which 
A 7^ 0. If we restrict our attention to initial data for which A < everywhere, then the evolution of the Ai 
would be determined and the Hamiltonian would be positive. We have 



H=-j£x (^-4. -I- {d^A^f ~ {d^A^f ) (102) 
= \j d^'x + F^,Fo, - 2{d,A^)F,^ ] (103) 

= ^ y d'a; Q^^f, + F^,Fo, + 2Aod,F,o ] (104) 



2 J d-'x ^-Ft^ + Fo^Fo, ~ 4AA^ j , (105) 

which is manifestly positive when A < 0. However, it is not clear why we should be restricted to this form 
of initial data, nor whether even this restriction is enough to ensure stability beyond perturbation theory. 

The status of this model in both the spacelike and timelike cases remains unclear. However, there are 
indications of further problems. For the spacelike case, Peloso tt. al. find a linear instability for perturbations 
with wave numbers on the order of the Hubble parameter in an exponentially expanding cosmology [34, 3^. 
For the timelike case, Seifert found a gravitational instability in the presence of a spherically symmetric 
source EBi. 



B. Scalar action 



The equation of motion for the scalar Lagrangian with a fixed-norm constraint is 

d-'d^A^ = 2XA''. 



(106) 
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Using the fixed-norm constraint {Afj^Ai^ = ), we can solve for the Lagrange multiplier field, 

A = T^A.5''5^A^ (107) 

In contrast with the Maxwell theory, in the scalar theory it is the timelike case for which we can demonstrate 
obstacles to smooth evolution, while the spacelike case is less clear. (The Hamiltonian is bounded below, 
but there are no perturbative instabilities or known obstacles to smooth evolution.) 
When the vector field is timelike, we have four constraint equations in the scalar case, 

Al - AiAi = and di{d^Ai^) = 2XAi. (108) 

Suppose we give homogeneous initial conditions such that ^o(io) > ti- Align axes such that, 

AM(^o) = (Ao(to),0,0,A3(to)), (109) 

where A^it^Y = Aoit^Y — m^. Note that, since ^3(^0) 7^ 0, we have that A = from the v = 3 equation of 
motion. The v = equation of motion therefore gives, 

(no, 

We see that the timelike component of the vector field has the time-evolution, 

Ao{t) = Ao{to) + C{t-to). (Ill) 

For generic homogeneous initial conditions, C 7^ 0. In this case, Aq will not have a smooth time evolution 
since Aq will saturate the fixed-norm constraint, and beyond this point Aq will continue to decrease in 
magnitude. To satisfy the fixed-norm constraint, the spatial components of the vector field Ai would need 
to be imaginary, which is unacceptable since is a real-valued vector field. This problem never occurs for 
the spacelike case since there always exist real values of Ai that satisfy the c;oiistraint for any Aq. 

Again, it is possible that this obstruction to a well-defined evolution will be regulated by terms of higher 
order in the effective field theory. The time derivative of A3 is 

^''^s = f° ., d^Ao. (112) 

As AqAq approaches mP, with finite derivatives of the derivative of the spatial component A3 becomes 
unbounded. If higher-order terms in the effective action have time derivatives of the components Ai, these 
terms could become relevant to the vector field's dynamical evolution, indicating that we have left the realm 
of validity of the low-energy effective field theory we axe considering. 

Whether or not a theory with a scalar kinetic term and fixed expectation value is viable remains uncertain. 



VII. CONCLUSIONS 



In this paper, we addressed the issue of stability in theories in which Lorentz invariance is spontaneously 
broken by a dynamical fixed-norm vector field with an action 

d^x i--p,F^,F''^ - Md^A^^ - i3^—^{d^Ap){d,AP) + XiA^^A^ ± m^) j , (113) 

where A is a Lagrange multiplier that strictly enforces the fixed-norm constraint. In the spirit of effective 
field theory, we limited our attention to only kinetic terms that are quadratic in derivatives, and took care 
to ensure that our discussion applies to regimes in which an effective field theory expansion is valid. 

We examined the boundedness of the Hamiltonian of the theory and showed that, for generic choices of 
kinetic term, the Hamiltonian is unbounded from below. Thus for a generic kinetic term, we have shown that 
a constant fixed-norm background is not the true vacuum of the theory. The only exception is the timelike 
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sigma-model Lagrangian (/3i — P^,, f34 — and A^^A^ = — m^), in which case the Hamiltonian is positive- 
definite, ensuring stabihty. However, if the vector field instead acquires its vacuum expectation value by 
minimizing a smooth potential, we demonstrated (as was done previously in |14j ) that the theory is plagued 
by the existence of a tachyonic ghost, and the Hamiltonian is unbounded from below. The timelike fixed- 
norm sigma-model theory nevertheless serves as a viable starting point for phenomenological investigations 
of Lorentz invariance; we explore some of this phenomenology in a separate paper [22j . 

We next examined the dispersion relations and energies of first-order perturbations about constant back- 
ground configurations. We showed that, in addition to the sigma-model case, there are only two other choices 
of kinetic term for which perturbations have non-negative energies and do not grow exponentially in any 
frame: the Maxwell — P4 — 0) and scalar (/3i — (3^ — 0) Lagrangians. In either case, the theory has fewer 
than three propagating degrees of freedom at the linear level, as some of the vector components in the action 
lack time derivatives and act as additional Lagrange multipliers. A subset of the phase space for the Maxwell 
theory with a timelike aether field is well-defined and stable, but is identical to ordinary electromagnetism. 
For the Maxwell theory with a spacelike aether field, or the scalar theory with a timelike field, we can find 
explicit obstructions to smooth time evolution. It remains unclear whether the timelike Maxwell theory or 
the spacelike scalar theory can exhibit true violation of Lorentz invariance while remaining well-behaved. 
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APPENDIX A: SOLUTIONS TO THE LINEARIZED EQUATIONS OF MOTION 

We start by finding the solution to the equations of motion, linearized about a timelike, fixed-norm 
background, Then, showing less details, we find the solutions to the equations of motion linearized 



about a spacelike background. Finally, we put the solutions in both cases into the compact form of (A26) 



(A28|. Our results agree with the solutions for Goldstone modes found in 



The equations of motion for a timelike (-I-) or spacelike (— ) vector field are (16), 



Q,, = ) {(iidpdPA- + - ii^)d^dpAf + fi^G^) = 0, (Al) 



where is defined in (14 1 and A^Q^ = identically. 

Timelike background. Consider perturbations about an arbitrary, constant (in space and time) timelike 
background = that satisfies the constraint: Af^A'^ = —ni^. Define perturbations by = + 
Then, to first order in these perturbations, A^Q^ = identically, and r]^" A^SAi, = by the constraint. We 
can define a basis set of four Lorentz 4- vectors , with components 

nl = Ajm, 71; ; ze{l,2,3}, (A2) 

such that 

Ty^-^n^^ = (A3) 

The independent perturbations are (5a" = •q^'^n'^SAi, for a = 1,2,3. {5aP is zero at first order in pertur- 
bations due to the constraint.) It is then cleax that there are three independent equations of motion at first 
order in pertubations (assuming the constraint) for the three independent perturbations, 

5Q' = nl {iSidpdPSA" + {(3, - Pi)d''dpdA'' + p^nln^d^'dPdA") = 0, (A4) 

where i S {1, 2, 3}. We look for plane wave solutions for the 6 A: 

5Ap= f d''kqpik)e'''"^\ (A5) 
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Since ij^'^n^SAi, — 0, at first order, 

1t^ = '^jK ^^^'^^ je {1,2,3}. (A6) 

The equations of motion become the algebraic equations: 

= {PikpkPnln^'' + (/3, - /3iXfc''n^fc'^ + /34n^nOfc^A:^n>^''') 9 (A7) 

= {PikpkPS'^ + {(5, - (3i)nlk''n'^k'' + l3inlnlk^'kP5'^) Cj (A8) 

= M'^Cj. (A9) 

The three independent solutions to these equations are given by setting an eigenvalue of the matrix M 
to zero and setting Ci to the corresponding eigenvector. Setting an eigenvalue of M equal to zero gives a 
dispersion relation, 

l3ikpkP + l3i{nlkPf = (AlO) 
with two linearly independent eigenvectors, 

(62)* = t2ijn'pkP ; (63)^ = Cs.jnlkP. (All) 
The second eigenvalue of M gives the dispersion relation, 

P.kpkP + - /3i + /34)Kfc'^)' - 0, (A12) 

with corresponding eigenvector, 

c, = n'pkP. (A13) 

Spacelike background. The first order linearized equations of motion about a spacelike background are: 

SQ" = nl {pidpdPSA'' + {P, - (3i)d''dpSAP + PinlnldPdPSA") = (A14) 

where a € {0, 1, 2} and where, similarly to the timelike case, we have defined the set of four Lorentz 4- vectors, 
n", to be 

nl = Ap/m and n^;ae{0,l,2} (A15) 

such that 

^fripu'^ = V""^- (A16) 

The independent perturbations are da°' = -qP"" A,y for a = 0, 1, 2. {5a^ is zero at first order in perturbations 
due to the constraint.) 

Again we look for plane wave solutions of the form in (A5l. But now, since rjP'^n^SAi, — 0, at first order, 

Qp = Can'^, where a e {0, 1, 2}. (A17) 
The equations of motion become the algebraic equations: 

= {pikpkPny'"' + iP, - P^)nlVn\kP + P^nl'nlkPkPnln^'') Cf, (A18) 
= {PikpkPi^'^'' + {13, - l3i)nlk''nlkP + P^nlnlkPkPr^'''') Cb (A19) 
= M°'''cb. a, &e {0,1,2} (A20) 
Two independent solutions correspond to the dispersion relation (a e {0, 1, 2}) 

(3ikpkP + (3i{nlkPy ^ , (A21) 

with corresponding eigenmodes 

{ei)a = eaib3nlkP ; (e2)a = eab23<fc''. (A22) 
The third solution corresponds to the dispersion relation 

P,kpkP - iP, - /3i - P4){nlkPf = , (A23) 

with corresponding eigenmode 

Ca = VabulkP. (A24) 
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General expression. We can express the solutions in the timelike and spacelike cases in a compact form 
by using the orthonormality of the n", (A3 1, along with (A2|, (A15I, and the fact that/ 



Then plugging (A6l and (A17l into (A5) yields the solutions, 
where either, 

qJk) = ia^kP — e„^„o- and PikpkP + (5^ — — 
m \ m 

where are real-valued constants or, 



and a^A^ = 0, 



(A25) 



(A26) 



(A27) 



and 



p.kpkP ± (13, - f3i ± Pi) 



ApkP 



(A28) 



where a is a real-valued constant. The reality of the a's follows from the condition, qp{k) = q*^{^k), that 
holds if and only if 5Ap in ( |A5p is real. In (A28l, the sign corresponds to the timelike background and 
the "— " sign to a spacelike background. 
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